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Abstract. We performed a deep infrared imaging survey of 63 embedded young stellar objects (YSOs) located in 
the Taurus and Ophiuchus clouds to search for companions. The sample includes Class I and flat infrared spectrum 
protostellar objects. We find 17 companions physically bound to 15 YSOs with angular separations in the range 
0.8-10" (110-1400 AU) and derive a companion star fraction of 23±9% and 29±7% for embedded YSOs in Taurus 
^ . and Ophiuchus, respectively, about twice as large as that found among G dwarfs in the solar neighborhood. 

' Therefore, binary and multiple protostellar systems are a very frequent outcome of the fragmentation of prestellar 

QQ , cores. In spite of different properties of the clouds and especially of the prestellar cores, the fraction of wide 

■ companions, 27±6 % for the combined sample, is identical in the two star-forming regions. This suggests that the 
' frequency and properties of wide multiple protostellar systems are not very sensitive to specific initial conditions. 
, Comparing the companion star fraction of the youngest YSOs still surrounded by extended envelopes to that 

JnI ■ of more evolved YSOs, we find evidence for a possible evolution of the fraction of wide multiple systems, which 

seems to decrease by a factor of about 2 on a timescale of ~10^ yr. For the first time, it is possible to confront 
. the result of a multiplicity survey of a nearly complete population of embedded YSOs at an age of ~10^yr to 

CL(' numerical simulations of molecular cloud collapse which, after a few free fall times, reach this evolutionary stage. 

JL I Somewhat contrary to model predictions, we do not find evidence for a sub-clustering of embedded sources at this 

■ stage on a scale of a few 100 AU that could be related to the formation of small-N protostellar clusters. Possible 
' interpretations for this discrepancy are discussed. 

> 
X 

^ 1. Introduction is currently under debate. Numerical simulations suggest 

■ - - that the dynamical fragmentation process is strongly de- 

Giant molecular clouds produce stars in a multi-step pro- pendent on the initial conditions, such as cloud rotation, 
cess that needs to be understood m order to account temperature or magnetic field strength (e.g., Bonnell 

for the general properties of stars, such as the steUar ^gg^. j^^^j^^^ ^ ^^^^^-^ ^994. g^^^ 2002). With 

mass function or the properties of muhiple systems m the ^^^^^^1 star-forming regions located within 500 pc of the 
Galaxy. First, the molecular cloud fragments into mdi- ^^^^ ^.^^^^ ^^^^^ relatively quiescent zones such as the 

vidual cores that subsequently collapse to form individual Taurus- Auriga cloud to the dense cluster-forming Orion 
stars. In many cases, however, those cores fragment a sec- Nebula, a number of observational tests can be performed 
ond time during protostellar collapse (m the following, we ^^^^^-^ environment-induced effects, 
refer to this phenomenon as "core dynamical fragmenta- 

tion"), which leads to the formation of binary and higher ^he study of the collapse and core fragmenta- 
order multiple systems. The extent to which this pro- ^^"'^ Processes usually focuses on two main approaches. 

J J • -J-- 1 1, • 1 J -J-- r j-i 1 J Investigation of prestellar cores within molecular clouds 

cess depends upon initial physical conditions 01 the cloud ... 

provides insight on the initial conditions and the very early 

W offprint ^sts to: G. Duchene, e-mafl: phases of the collapse while the analysis of populations of 

gaspard.duchene@obs.ujf-grenoble.fr pre-main sequence (TTauri) objects reveal their outcome. 

* Based on observations obtained at the Canada-France- Due to the Umited spatial resolution of the single-dish mil- 
Hawaii Telescope and at the European Southern Observatory, limeter images from which prestellar core properties are 
La Silla, Chile. derived, it is usually not possible to probe spatial scales of 
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a few tens or hundreds AU, typical of most multiple sys- 
tems. On the other hand, TTauri stars have not necessar- 
ily preserved their original properties after the dissipation 
of their collapsing envelope. Studying more deeply embed- 
ded young stellar objects (YSOs), which are assumed to 
be much younger than T Tauri stars, at high spatial res- 
olution (< 1") offers an opportunity to sample a crucial 
intermediate stage between prestellar cores and pre-main 
sequence stars. 

Recent studies have revealed ranges in the properties 
of prestellar cores as well as in the young populations re- 
sulting from the star formation process, with significant 
differences between star-forming regions. We summarize 
here some of the most notable ones, with emphasis on 
the Taurus and Ophiuchus molecular clouds which con- 
stitute the focus of the present study. First of all, the 
Taurus molecular cloud forms only a few tens of objects in 
~ 1 pc-radius groups (Gomez et al. 1993) and is frequently 
referred to as the prototype of "isolated" star formation. 
It is believed that most stars form in much richer, cluster- 
like structures with several hundred objects within a 1 pc- 
radius (e.g.. Carpenter 2000; Adams & Myers 2001). While 
the Orion Trapezium cluster is a well known example of 
this "clustered" mode of star formation, the Ophiuchus 
molecular cloud also belongs to this category, as it host 
several hundred objects with an average spatial density 
on the order of 10^ stars/pc^ (Bontemps et al. 2001; Allen 
et al. 2002). Taurus and Ophiuchus are therefore excel- 
lent candidates to investigate whether the star formation 
process is sensitive to environmental conditions. 

Millimeter continuum mapping of both star-forming 
regions have revealed clear differences in the properties 
of their prestellar cores. While cores in Taurus appear to 
have a radius as large as 10000-20000 AU and smoothly 
merge into the ambient molecular cloud (Motte & Andre 
2001), those found in Ophiuchus exhibit a sharp bound- 
ary at a radius on the order of 5000 AU (Motte, Andre & 
Neri 1998), possibly due to external perturbations in this 
denser environment. The derived envelope masses, bolo- 
metric luminosities and density profiles of prestellar cores 
also differ between the two molecular clouds. Likewise, 
at a slightly later stage, the circumstellar envelopes of 
Ophiuchus embedded protostars are smaller than those 
of their Taurus counterparts (Motte et al. 1998). This led 
Motte & Andre (2001) to state that "the beginning of pro- 
tostellar evolution is suggested to be more violent (...) in 
protoclusters compared to regions of distributed star for- 
mation like Taurus" . In other words, the initial conditions 
for core fragmentation appear to be significantly different 
in the Taurus and Ophiuchus molecular clouds. 

Multiple systems are useful tracers of the core frag- 
mentation process. Systematic surveys have revealed that 
young stellar populations in dense stellar clusters (Orion 
Trapezium, IC 348, Pleiades) exhibit multiplicity rates 
similar to field stars whereas populations of T Tauri stars 
located in more distributed associations can possess up to 
twice as many companions (e.g., Duchene 1999). However, 
it must be noted that the frequency of visual companions 



to T Tauri stars in the Taurus and Ophiuchus clouds are 
identical, both exhibiting a strong excess of multiple sys- 
tems, i.e., that the latter behaves in the same way as "loose 
associations" in this respect although it is an example of 
clustered star formation. 

The lower multiplicity rate in dense clusters does not 
necessarily imply that the core fragmentation process it- 
self is environment-dependent. For instance, Kroupa, Petr 
& McCaughrean (1999) have performed N-body simula- 
tions to follow the dynamical evolution of young stellar 
clusters and showed that the observed properties of mul- 
tiple systems in these populations are consistent with an 
initial set of multiple systems identical to loose T Tauri 
associations. In clusters, direct dynamical encounters be- 
tween systems are able to disrupt a fraction of the ini- 
tial (wide) systems over timescales as short a 1 Myr (e.g., 
Kroupa 1995a), thus rapidly bringing the multiplicity rate 
down to the observed value. 

In order to probe the outcome of the fragmentation 
process before it is dynamically modified, at least in clus- 
tered environments, it is therefore necessary to study 
younger YSOs than T Tauri stars. In this study, we present 
a systematic search for close companions among Class I 
sources and slightly more evolved "flat-spectrum" objects 
in Taurus-Auriga and Ophiuchus. These objects exhibit a 
flat or rising spectral energy distribution from the near to 
the mid infrared and are believed to still be heavily em- 
bedded within a collapsing envelope. On a statistical basis, 
their estimated age is of order of 10^ yrs (e.g. Kenyon & 
Hartmann 1995). Because they are seen through a high 
column density of dust, these objects are almost unde- 
tected in the visible and deep near infrared imaging is 
required to look for close companions. Recently, Haisch 
et al. (2004) reported the results of an imaging survey 
of samples of embedded YSOs located in several nearby 
molecular clouds and derived a multiplicity rate over a sep- 
aration range 300-2000 AU similar to that of more evolved 
T Tauri stars but they did not conclude regarding possi- 
ble differences between various star-forming regions. Our 
survey differs from theirs in two respects. Firstly, we focus 
on homogeneous and nearly complete samples of embed- 
ded YSOs in two star-forming regions only, Taurus and 
Ophiuchus, in order to be able to contrast at a significant 
statistical level the multiplicity rate resulting from the iso- 
lated and clustered star formation modes. Second, the 0.6" 
angular resolution of the present survey allows us to ex- 
plore closer systems down to a separation of 110 AU at the 
distance of these two clouds (140 pc, Bertout, Robichon 
& Arenou 1999; Bontemps et al. 2001). 

This paper is organized as follows: in §[2 we present 
our samples and deep infrared observations obtained at 
CFHT and ESO/NTT. Candidate visual companions de- 
tected on the images are presented in §|31and physical com- 
panions to embedded YSOs are identified on statistical 
grounds, thus allowing us to establish the multiplicity fre- 
quency of embedded YSOs in the Taurus and Ophiuchus 
clouds. In §01 we discuss the results and their implica- 
tions regarding the fragmentation of prestellar cores in the 
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Table 1. Sample of protostars surveyed in Taurus. The infrared photometry and the spectral indices are taken from 
Kenyon & Hartmann (1995), except for IRAS 04191+1523, whose JHK photometry is from the 2MASS database and 
classification from Motte & Andre (2001). The bolometric luminosities are from Motte & Andre (2001). IRAS did not 
resolve the systems IRAS 04108+2803 AB and IRAS 04181+2654 AB and we assigned the same spectral index, based 
on unresolved photometry, to both components and split evenly the bolometric luminosity of the later system. The 
classification of each source is based on the infrared spectral indices (see text); "I" stands for Class I sources and "FS" 
for flat spectrum sources. The last three columns indicate whether the source i) has at least one visual companion 
within 10", ii) is surrounded by an extended nebulosity in our images (see also Park & Kenyon 2002), and iii) has an 
extended (> 1000 AU) envelope as seen in millimeter continuum mapping (from Motte & Andre 2001). 



Object 


J 


H 


K 






Class 


Lboi{LQ) 


Comp? 


Opt neb? 


Mm env? 


IRAS 04016+2610 


13.62 


11.67 


9.33 


1.06 


1.02 


I 


3.7 


N 


Y 


N 


IRAS 04108+2803 B 


16.25 


13.25 


11.12 


0.50 


0.62 


FS 


0.6 


N 


N 


N 


IRAS 04113+2758 


11.29 


9.15 


7.79 


-0.01 


0.12 


FS 


> 1.6 


Y 


N 


Y 


IRAS 04158+2805 


13.15 


11.92 


10.80 


1.47 


0.72 


I 


> 0.4 


N 


Y 


N 


IRAS 04169+2702 


16.87 


13.83 


11.22 


0.89 


1.10 


I 


0.8 


N 


Y 


Y 


IRAS 04181+2655 


15.46 


12.74 


10.54 


0.60 


0.55 


I 


0.4 


N 


N 


N 


IRAS 04181+2654 A 


15.90 


12.74 


10.58 


-0.00 


0.10 


FS 


0.25 


N 


N 


Y 


IRAS 04181+2654 B 


18.11 


13.64 


10.81 


-0.00 


0.10 


FS 


0.25 


N 


N 


Y 


IRAS 04191+1523 


16.74 


14.43 


12.26 






I 


0.5 


Y 


Y 


Y 


IRAS 04239+2436 


15.61 


12.96 


10.58 


1.27 


1.15 


I 


1.3 


N 


N 


Y 


IRAS 04248+2612 


12.76 


11.39 


10.65 


-0.09 


0.48 


FS 


0.4 


Y 


Y 


Y 


IRAS 04260+2642 


14.64 


12.96 


11.63 


0.16 


0.25 


FS 


0.1 


N 


N 


N 


IRAS 04263+2426 


10.84 


9.06 


7.80 


1.02 


0.74 


I 


7.0 


Y 


Y 


Y 


IRAS 04264+2433 


13.63 


12.07 


11.09 


0.67 


0.98 


I 


0.5 


N 


N 


N 


IRAS 04295+2251 


13.62 


11.11 


9.55 


0.11 


0.21 


FS 


0.6 


N 


N 


N 


IRAS 04302+2247 


14.20 


12.35 


11.07 




0.15 


FS 


0.3 


N 


Y 


N 


IRAS 04325+2402 




13.02 


11.03 




0.79 


I 


0.9 


Y 


Y 


Y 


IRAS 04361+2547 


16.07 


12.88 


10.55 


1.27 


1.55 


I 


3.7 


N 


Y 


Y 


IRAS 04365+2535 


17.28 


13.39 


10.62 


1.08 


1.27 


I 


2.4 


N 


Y 


Y 


IRAS 04381+2540 


16.21 


14.32 


12.00 


1.20 


1.31 


I 


0.7 


N 


Y 


Y 


IRAS 04385+2550 


12.61 


10.72 


9.65 


0.14 


0.18 


FS 


> 0.4 


N 


N 


N 


IRAS 04489+3042 


13.74 


11.58 


10.11 


0.18 


0.31 


FS 


0.3 


N 


N 


N 



two star- forming regions, and compare the multiplicity fre- 
quency we derive for embedded Class I and flat spectrum 
sources to that of both younger Class sources and older 
Class II-III T Tauri stars in the same clouds in order to in- 
vestigate the dynamical evolution of protostellar systems 
on a timescale of a few Myr. Finally, we summarize our 
results in §[S1 

2. Sample selection and observations 

We have built our samples from the object lists from 
Motte & Andre (2001) in Taurus and Bontemps et al. 
(2001) in Ophiuchus. The list of our targets is presented in 
Tables^ and 121 together with their near-infrared photom- 
etry and near- to mid-infrared spectral index. The latter 
is used to distinguish between Class I, flat-spectrum and 
T Tauri (Class II) sources, where we adopted am — 0.55 
and am — —0.05 thresholds between these categories (see 
Bontemps et al. 2001). Note that because of the large line- 
of-sight extinctions towards some of our targets, the use of 
the observed infrared spectral index is not completely se- 
cure. While there might be a few embedded T Tauri stars 
in our sample, especially among the objects classified as 
flat spectrum, we consistently use the observed spectral 
index to estimate each object's nature since we cannot es- 



timate the extinction to each individual source. Overall, 
we have observed all known sources in both areas with the 
exception of two Class I sources and one flat-spectrum ob- 
ject in Taurus, all of which are in the crowded vicinity of 
HL Tau, to the south of the molecular cloud. Our sample 
therefore includes 22 protostellar sources in Taurus and 
41 in Ophiuchus, in the magnitude range 7.5 < K < 14.5 
with a median of K ^ 11 mag. A near-infrared color-color 
diagram of all our targets is presented in Figure ^ which 
shows that most sources have redder colors than T Tauri 
stars, even if extinction is taken into account. 

We conducted most of our observations with the wide- 
field infrared camera CFHT-IR (Starr et al. 2000) at 
Canada- France-Hawaii Telescope, a 1024x1024 detector 
which offers a pixel scale of 0'.' 211 /pixel and a total field-of- 
view of 3.'6. The Ophiuchus images were obtained on June, 
11th and 12th 2001, while the Taurus- Auriga dataset was 
obtained on October, 29th and 30th 2001. Each target was 
systematically surveyed in the K band, where these proto- 
stars are brighter. In some multiple systems, we obtained 
follow-up H and J images to compare the near-infrared 
colors of the primary to its candidate companions. For 
each field, a sequence of 4 images was obtained, each of 
them being cubes of 4 exposures. Individual integration 
times ranged from 0.1s to 40 s at if and up to Imin at 
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Table 2. Sample of protostars surveyed in Ophiuchus. ISO numbers, infrared spectral indices and bolometric lumi- 
nosities are from Bontemps et al. (2001), while the near-infrared magnitudes are from Barsony et al. (1997), Greene 
et al. (1994) and the 2MASS database (2MASS uses a Ks filter rather than a K filter; the difference in magnitude 
for an object can reach a few tenths of a magnitude). The source marked with a f symbol has an infrared spectral 
index typical of a normal Class II source but Ressler & Barsony (2001) have shown that this system contains a Class I 
source and it is classified as such here. The last three columns have the same meaning as in Tabled The existence of 
an envelope is taken either from Andre & Montmerle (1994) and Motte et al. (1998) or from our reexamination of the 
mosaic presented in the latter. (Note that this classification is more uncertain than in Taurus due to the smaller size 
of protostellar envelopes in Ophiuchus. Particularly uncertain cases are flagged with a "?".) 



ISO# 


Object 


J 


H 


K 


•2-14, 
(^IR 


Class 




Comp? 


Opt neb? 


Mm env? 


21 


CRBR 12 


>17.00 


15.58 


12.04 


0.91 


I 


0.42 


N 


N 


N 


26 


CRBR 15 


16.35 


13.89 


11.94 


0.01 


FS 


0.083 


N 


N 


N 


29 


GSS30, GY6 


13.89 


10.83 


8.32 


1.20 


I 


21. 


N 


Y 


Y? 


31 


LFAMl 


>17.00 


>15.50 


13.59 


1.08 


I 


0.13 


Y 


N 


Y 


33 


GYll 


16.52 


15.37 


14.15 


0.31 


FS 


0.011 


Y 


Y 


N 


46 


VSSG27, GY51 


16.69 


13.46 


10.72 


0.17 


FS 


0.41 


Y 


N 


N 


54 


GY91, CRBR 42 


>17.00 


16.40 


12.51 


0.70 


I 


0.17 


N 


N 


Y? 


65 


WL12, GYlll 


16.86 


13.07 


10.18 


1.04 


I 


2.6 


N 


Y 


Y? 


70 


WL2, GY128 


>17.00 


14.05 


10.99 


0.05 


FS 


1.4 


Y 


N 


N 


75 


GY144 


>17.00 


15.75 


13.46 


0.20 


FS 


0.053 


N 


N 


N 


77 


GY152 


>17.00 


>15.75 


13.64 


0.05 


FS 


0.052 


N 


N 


N? 


85 


CRBR 51 


>17.00 


>15.75 


14.00 


0.03 


FS 


0.035 


N 


N 


N? 


99 


LFAM 26, GY 197 


>17.00 


>17.90 


14.63 


1.25 


I 


0.064 


N 


Y 


Y 


103 


WL17, GY205 


16.90 


13.57 


10.28 


0.42 


FS 


0.88 


N 


N 


Y? 


108 


EL 29, GY214 


17.21 


12.01 


7.54 


0.98 


I 


26. 


N 


Y 


Y 


112 


GY224 


>17.00 


13.69 


10.79 


0.34 


FS 


1.7 


N 


N 


N? 


118 


IRS 33, GY236 


>17.00 


15.27 


12.23 


0.32 


FS 


0.43 


N 


N 


N? 


119 


IRS 35, GY238 


>17.00 


16.25 


12.73 


0.30 


FS 


0.45 


N 


N 


N 


121 


WL20, GY240 f 


13.45 


10.78 


9.21 


-0.07 


I 


1.5 


Y 


N 


Y? 


124 


IRS 37, GY244 


>17.00 


13.90 


10.95 


0.35 


FS 


1.5 


Y 


Y 


N 


127 


GY245 


19.00 


15.43 


11.98 


0.17 


FS 


0.13 


N 


N 


N 


129 


WL3, GY249 


>17.00 


14.49 


11.20 


0.23 


FS 


1.6 


N 


Y 


N? 


132 


IRS 42, GY252 


15.21 


11.31 


8.41 


0.08 


FS 


5.6 


N 


N 


N? 


134 


WL6, GY254 


>17.00 


14.39 


10.04 


0.59 


I 


1.7 


N 


N 


N 


137 


CRBR 85 


>17.00 


17.13 


13.21 


1.48 


I 


0.36 


N 


N 


Y 


139 


GY 260 


>17.00 


15.71 


12.54 


-0.03 


FS 


0.30 


N 


N 


N 


141 


IRS 43, GY265 


>17.00 


13.17 


9.46 


0.98 


I 


6.7 


Y 


Y 


Y 


143 


IRS 44, GY269 


>17.00 


13.09 


9.65 


1.57 


I 


8.7 


N 


Y 


Y 


145 


IRS 46, GY274 


>17.00 


14.65 


11.46 


0.94 


I 


0.62 


N 


N 


Y? 


147 


IRS 47, GY279 


15.44 


11.64 


8.95 


0.17 


FS 


3.7 


N 


Y 


N 


159 


IRS 48, GY304 


10.53 


8.65 


7.42 


0.18 


FS 


7.4 


N 


N 


N 


161 


GY301 


>17.00 


14.77 


11.30 


0.12 


FS 


1.8 


N 


N 


N? 


165 


GY312 


16.23 


13.75 


11.93 


0.03 


FS 


0.091 


N 


N 


N 


167 


IRS 51, GY315 


17.12 


12.42 


8.93 


-0.04 


FS 


1.1 


Y 


N 


Y? 


170 


B 162741-244645 


>17.00 


15.71 


13.54 


0.51 


FS 


0.065 


N 


N 


N 


175 


GY344 


16.81 


14.05 


11.84 


0.02 


FS 


0.10 


N 


Y 


N 


182 


IRS 54, GY378 


16.63 


13.50 


10.87 


1.76 


I 


6.6 


Y 


Y 


Y? 


200 


ISO 1631434-245524 


14.31 


12.00 


10.43 


0.22 


FS 


1.3 


N 


N 


N 


203 


ISO 1631520-245536 


>18.45 


14.86 


12.73 


1.19 


I 


1.0 


N 


Y 


N? 


209 


IRS 67 


>17.30 


13.31 


10.43 


0.75 


I 


1.5 


Y 


Y 


Y 


210 


ISO 1632021-245616 


15.81 


14.08 


13.02 


0.09 


FS 


0.094 


Y 


N 


N 



H and J. In some Ophiuchus fields, we obtained two sets 
of images at the same location, with a short and a long 
integration time in order to observe with similar contrasts 
both bright and faint targets in crowded fields. The me- 
dian value and standard deviation of the full width at half 
maximum (FWHM) measured on a sample of single stars 



throughout the runs are 0'.'60±0'.'07 and 0'.'65±G'.'1G for the 
Ophiuchus and Taurus datasets respectively. 

On March, 10th 2000, we obtained additional images 
of the L 1689 area, a core located to the south of the main 
cores of the Ophiuchus molecular cloud, with ESO's New 
Technology Telescope at la Silla Observatory. We used 
SOFI, which has a 1024x1024 detector, a O'.'288/pixel 
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Fig. 1. Near-infrared color-color magnitude for all our 
targets in Taurus {filled triangles) and Ophiuchus (filled 
squares). Objects not detected by 2MASS in the J band 
are indicated with arrows. The companions detected at 
all three wavelengths in our survey are indicated by dia- 
monds, the larger ones representing those that are likely 
to be bound to one of our targets. The solid curve indi- 
cates the locus of main sequence sources and the dashed 
line represents the dereddened locus of all TTauri stars 
in the Taurus region (Meyer 1996). The dotted lines de- 
lineate the reddening band of both loci for Ay — 25 mag 
and the extinction law from Rieke & Lebofsky (1985). 

scale and 4.'9 field-of-view centered on the known T Tauri 
star L1689-IRS5. With all three J, H and K filters, we 
obtained a sequence of 10x1.18s exposures using a ran- 
dom jittering pattern in a 20"x20" box. The final images 
contain two Class I sources and two flat spectrum objects. 
The image quality was comparable to that of our CFHT 
images, with a 0'.'55 FWHM, which enables us to combine 
this image to our survey. 

All images of a given object with the same integration 
time were median-combined to obtain a sky estimate, that 
was then subtracted from each individual image. The im- 
ages were then flat-fielded, aligned based on the location 
of the brightest unsaturated source in the field, and me- 
dian combined to produce the final images used in this 
study. The astrometry of the multiple systems was esti- 
mated using the known pixel scale and orientation of the 
detectors; uncertainties are on the order of 0.5% for the 
plate scale and 1° for the orientation. 

3. Results 

The objective of our study is to identify physical compan- 
ions to embedded protostars in Taurus and Ophiuchus. We 
first identify all visual companions in the 0.8-10" (110- 
1400 AU) separation range fS I3.1|l and compare our re- 
sults to those of previous studies of protostars in both 
areas fS I3.2(l . We then estimate which of these compan- 
ions are likely to be merely chance projections of back- 



ground/foreground objects and finally derive the resulting 
multiplicity frequency of embedded YSOs in Taurus and 
Ophiuchus fSrOll. 

3.1. Visual companions to Class I and flat spectrum 
sources 

Although embedded sources are located in dense molecu- 
lar clouds, imaging studies as deep as the one presented 
here can reveal a large number of faint objects that are 
mostly unrelated to the star-forming region. We there- 
fore limited the search for companions to within a 10" 
(1400 AU) radius around the surveyed embedded YSOs, 
a cutoff which excludes the presence of many background 
stars projected by chance in the vicinity of out targets (see 
Using this search radius also allows to compare di- 
rectly our results to multiplicity surveys among the more 
evolved T Tauri stars, which have been conducted over the 
same separation range. Furthermore, if the distribution of 
orbital periods of protostellar multiple systems does not 
evolve significantly with time (e.g., Kroupa & Burkert 
2001), we expect that this separation range will encom- 
pass a large fraction of all binary systems, as is the case 
for main sequence field stars: Duquennoy & Mayor (1991) 
have for instance shown that less than 10 % of solar-type 
field stars have a companion at separations wider than 
1400 AU. The lower limit on measurable binary separation 
is fixed by the sharpness of the images. Systems tighter 
than the image's FWHM would only be slightly elongated 
and could be misidentified as wind shaking or telescope 
tracking errors. In order to exclude binaries that are only 
marginally resolved, we do not consider systems tighter 
than 0'.'8 (110 AU). Tighter systems wiU be the focus of 
a higher resolution study of the same sources performed 
with adaptive optics (Duchene et al., in prep.). 

The 1400 AU upper cutoff adopted here is several times 
smaller than the size of individual prestellar cores in both 
star-forming regions (Motte et al. 1998, Motte & Andre 
2001). Therefore, two stars physically separated by such a 
"small" distance most likely formed through the fragmen- 
tation of a single core. However, it has been shown that 
the typical size of molecular cores differs in both areas, 
with cores in Taurus being ~ 3 times larger than those 
in Ophiuchus (see §P). In order to study the core frag- 
mentation process, it might be more relevant to search 
for binaries using a similar separation cutoff relative to 
the original core size. Therefore, we have also considered 
candidate companions up to 30" (4200 AU) from our YSO 
targets in Taurus, i.e., less than a third of the size of a typ- 
ical prestellar core in that region. In our analysis, we use 
the 10" cutoff to define our primary sample of compan- 
ions. The enlarged separation range probed in Taurus is 
only used for discussing core size-scaled binary properties 
in both star- forming regions (see §01 1). 

In the 0'.'8~-10" range, we found a possible companion 
to 5 of the 22 protostars we surveyed in Taurus, while in 
Ophiuchus, we found 14 candidate companions associated 
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Fig. 2. Greyscale images of the can- 
didate multiple systems detected in 
Opliiuchus. Each image is about 
20" (2800 AU) on a side; North is 
up and East to the left. The faintest 
companions are indicated with ar- 



to 11 targets in our sample of 41 protostars. The observed 
raw companion star fractions^ is 23 ±9% in Taurus and 
34 ±7% in Ophiuchus. In the 10-30" range around our 
Taurus targets, we find 19 additional possible compan- 
ions, although a majority of them are likely to be pro- 
jected background stars (see S I3.3|) . Images of the possible 
higher-order multiple systems with separations < 10" are 
presented in Figure 13.11 The astrometric and photomet- 
ric properties of all candidate companions are summarized 
in Table |31 The faintest two companions to IRS 54 were 
detected by combining a short integration image with a 
much longer exposure of the same field; they are therefore 
excluded from our discussion as we obtained such deep 
images for only a handful of objects in our sample. 

Our ability to detect a companion depends on its flux 
relative to its brighter primary: the closer the compan- 
ion, the brighter it needs to be for us to detect it. This 
is illustrated in Figure |31 where we have plotted all ob- 
served companions in our images. We have also included in 
this figure companions to well-exposed Ophiuchus T Tauri 
multiple systems that lie within our field of view (see 
Appendix^, in order to better define the lower enve- 
lope of all detected companions. It appears that binaries 
with Aif < 7-8 mag can be detected beyond 6" or so, 
while at 1", only binaries with AK < 3 mag are identi- 
fied. Separately, we have estimated our ability to detect 
faint point sources around bright sources by computing 
the standard deviation of pixel values in concentric annuli 
surrounding single sources. The resulting average 3 a limit 
is also plotted in Figure 13 along with the standard devia- 
tion (from observations of 16 independent point sources) 
about that curve. These curves match relatively well the 
lower envelope of our detected companions. Within the 
central 2-3", the detection limit mostly depends on the 
seeing but it does not appear to vary dramatically from a 
source to another. At larger distances, the detection limit 
is mostly driven by the background noise, which is roughly 
constant except if the protostar is surrounded by some ex- 
tended nebulosity. 

We conclude that we have detected all possible com- 
panions that lie above the solid curve in Figure 13 For ob- 
jects that have an extended nebulosity in our near-infrared 

^ Throughout this study we adopt the binomial statistics for 
estimating uncertainties associated with the companion star 
fraction, acsF = V ^comp(l — Ncomp/Nprim)/Nprim (whcre 
Nprim and Ncomp are the number of primaries and compan- 
ions, respectively), as a more appropriate method for this 
problem than the widely used Poisson statistics since the 
Ncomp/ Nprim ^ 1 Criteria is usually not met. 



images, our ability to detect a nearby companion can be 
somewhat impaired. As illustrated in Fi d3.2l nebulosities 
range from a faint extended halo (e.g., IRAS 04169-1-2702) 
to the two parallel nebulae associated with the presence of 
an edge-on, optically thick circumstellar disk that blocks 
the direct starlight (IRAS 04302-1-2247). However, only ex- 
tremely bright nebulosities in which no point-like object 
can be seen, such as IRAS 04302-1-2247, significantly affect 
our detection limit. Indeed, the fraction of visual compan- 
ions we detected around objects with an infrared nebu- 
losity (9 out of 24, 38 ±10%) is similar to that observed 
around point like sources (10 out of 39, 26 ±7%). We con- 
clude that we probably did not miss many companions 
due to the presence of extended nebulosities over the sep- 
aration range investigated in this study. 

3.2. Comparison to previous surveys 

In the Taurus- Auriga molecular cloud, this study is 
the first census of multiple systems among a nearly 
complete sample of embedded YSOs. We note that 
IRAS 04113-1-2758 is a resolved pair in 2MASS, while 
the binarity of IRAS 04263-^2426 (ee Haro6-10) has been 
discovered by Leinert & Haas (1989). The ~6" com- 
panion to IRAS 0419H-1523 has also been detected in 
the millimeter continuum with the IRAM Plateau de 
Bure Interferometer, confirming that both sources are 
young stellar objects with compact envelopes and/or disks 
(Motte 1998). Finally, our image reveal an almost point- 
like object located about 3" North of IRAS 04016-1-2610. 
Higher spatial resolution images of this system by Padgett 
et al. (1999) show an extended structure, most likely to be 
light scattered off the envelope or disk associated to this 
YSO. We therefore reject this "companion" from our anal- 
ysis. Of the 5 companions to Taurus embedded sources 
reported here, 2 are new detections (IRAS 04248-1-2612, 
IRAS 04325-f2402). 

Some of the widest pairs we identified in Ophiuchus 
were already known: LFAM 1 (whose companion is GY 12 
= ISO 34), IRS 43 (its companion is GY263), WL2 (from 
2MASS), WL20 (Ressler & Barsony 2001), GYll and 
IRS 37 (Allen et al. 2002). In their recent survey of em- 
bedded YSOs, Haisch et al. (2004) further identified the 
companions to IRS 54 and VSSG 27. The separations and 
position angle measurements for these systems agree well 
with their previously published values. We further report 
here the detection of companions to 2 targets that were 
also investigated by Haisch et al. (2002), namely IRS 51 
and IRS 67. Haisch et al. (2002) detected our proposed 



Duchene et al.: Multiple protostellar systems 



7 



Table 3. List of all companions detected within 10" (Ophiuchus) or 30" (Taurus) of our targets. Flux ratios are 
measured with respect to the brightest component in the K band image, except for three systems where the Class I 
source is the faintest component in the near-infrared (f symbol). The boldfaced entries represent systems that are 
likely to be physically bound (see 5 13. The two companions to IRS 54 listed in italics were detected by combining 
a short and a long exposure and are not included in our survey. Upper limits for flux ratios are at the 3 a level. 
IRAS 04264+2433 was saturated in all of our images, so only a lower limit for the flux ratio to its companion can be 
obtained. 
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companion to IRS 51 but their non-detection of the sec- 
ondary at 10 /im lead them to conclude that it was not 
real but rather a patch of scattered light in a possible 
outflow. The companion appears indistinguishable from a 
point source at the 0"6 resolution of our image of the sys- 
tem. We therefore believe that the companion is real and 
consider that its non-detection at 10 /xm results from the 
limited signal-to-noise ratio of the image of Haisch et al. 



(their Figure 1). The companion to IRS 67 was probably 
not detected by Haisch et al. because of its small separa- 
tion, 2"5, and large flux ratio, AK ^ 4.1 mag. Overall, 
of the 14 companions we report here for Ophiuchi em- 
bedded sources, 3 are new detections (IRS 51, IRS 67 and 
ISO 1632021-245616). 

Overall, our survey thus revealed 5 new companions 
in the 0"8-10" separation range, as well as the first in- 
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Fig. 3. Observed flux ratios of our targets (Taurus: trian- 
gles; Ophiuchus: squares) and of Class II and III systems 
in Ophiuchus (plus signs). The filled and empty symbols 
for our Ophiuchus targets discrimate between physically 
bound and projected systems (see i) l3.3|) . The sohd hue 
is our median Su detection limit, and the dashed lines 
represent the 1 a variation about that line (see text for 
details). 

frared detection of a companion only detected at mil- 
limeter wavelengths so far. This significantly adds to the 
statistics of multiplicity of embedded YSOs in both star- 
forming regions. Haisch et al. (2004) recently surveyed 6 
star- forming regions, with YSO samples ranging from 9 
to 30 objects per region, and derived an overall compan- 
ion star fraction of 12/76 (16 ±4%) over the separation 
range 300-1400 AU and down to a flux ratio AK=4.0 mag. 
Over the same separation range and detection limits, we 
derive a companion star fraction of 13/63 (20 ±5%) for 
embedded YSOs in Taurus and Ophiuchus. The 2 sur- 
veys yield very similar results, even though the former in- 
cludes smaller samples of YSOs distributed over a number 
of star-forming regions, while our survey deals with nearly 
complete YSO populations in only 2 star-forming regions. 
The similar results obtained from the two surveys thus 
suggests that the average companion star fraction for em- 
bedded sources does not vary much from one star-forming 
region to the other. Using high resolution radio continuum 
imaging, Reipurth et al. (1999, 2002, 2004) have compiled 
a list of 21 embedded protostars in several star- forming 
regions that includes 4 multiple systems with a separation 
in the range 110-1400 AU studied here, a multiplicity rate 
also very similar to the one we find here. 

3.3. Physically bound multiple systems 

The candidate companions we have identified could be ei- 
ther physically bound to our targets or unrelated objects 
(YSOs or field stars) that are projected within 10" (or 
30"). Since we only have K band photometry for most 
of our candidate companions, it usually is not possible to 



use color indices to determine their nature. The two visual 
companions to IRAS 04158-1-2805 are likely background 
stars, given their location in the near-infrared color-color 
diagram (see Figure^). In any case, the large and variable 
extinction across the molecular clouds would represent a 
major issue as a much redder companion could be inter- 
preted both as an extincted background star or a more 
embedded companion protostar. 

We therefore used a statistical approach to estimate 
the probability that the companions we found are phys- 
ically bound to their primary. Our method relies on the 
large field of view of our images, which allows us to esti- 
mate the local surface density of faint objects in the field. 
More specifically, we used the same method as described 
in Duchene et al. (2001), which we summarize here briefly. 

We first count the number of objects located in our 
field of view that are at least as bright as the candi- 
date companion, including the companion itself but not 
its primary. We convert this number into the average 
surface density of objects brighter than this threshold, 
I]{K < Kb). The median value of this quantity over 
all fields is about 2 x 10"'* arcsec"^. Assuming a ran- 
dom distribution of background stars across the flelds, we 
then estimated the number of companions one would ex- 
pect to lie as close to the primary as the companion or 
closer. This number is TV = Tr{p% ^ Pum) ^ ^ Kb), 
where pB is the binary separation and pum — 0'.'8 is the 
smallest separation probed in this survey. Finally, we es- 
timated the probability that a random uniform distribu- 
tion of stars across the field yields no companion within 
Pb, which we estimate using Poisson statistics (justified 
by the fact that N is almost always smaller than unity): 
Pbound — P{N = 0) = e~^. The resulting probabilities of 
the companion being physically associated to a YSO are 
given in the last column of Table 01 By construction, this 
method takes into account the projection of background 
and foreground objects, as well as of unrelated cloud mem- 
bers. 

Among the 0'.'8~10" companions detected in our sur- 
vey, 5 out of 14 in Ophiuchus and 3 out of 5 in Taurus 
have probabilities of being bound larger than 99 %. These 
are very likely physical systems but we have to decide on 
a threshold down to which we accept the companions as 
likely. On the other hand, most Taurus companions be- 
tween 10" and 30" have probabilities lower than 90 % of 
being bound. Those are likely to be projected unrelated 
stars. In the following, we decided to consider only systems 
that have probabilities of being bound larger than 95.5 %, 
i.e., for which a background status can be rejected at the 
2 (T level. This is probably a conservative estimate, as our 
method explicitly assumes that background stars are ran- 
domly located in the fields while our images show that 
the protostars are usually found in darker areas, with an 
increasing number of sources towards the edges of our im- 
ages (i.e., YSOs are located deeper in molecular clouds). 
This is however an image-to-image effect whose ampli- 
tude is difficult to estimate due to the relatively small 
number of stars per field. On the other hand, members of 
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Fig. 4. Example of nebulosities surrounding 
some of our targets; the upper row show some 
Taurus sources, the lower Ophiuchus targets. 
Each image is 25" on the side; North is up and 
East to the left. The contours are spaced by 
0.75 mag/arcsec^ from the peak flux down to 
approximately the 4-5 cr limit. 



young associations tend to cluster in relatively compact 
groups, which may produce apparent pairs of two YSOs 
that are in fact unbound. For instance, among the com- 
panions retained with this criterion, one (GY 12, compan- 
ion of LFAM 1) is a known member of the Ophiuchus star- 
forming region (see S I3.2|) that appears to be a Class III, 
weak-line T Tauri star. It is not clear whether two ob- 
jects that have such a different evolutionary status do form 
a physically bound binary system or whether this is the 
mere result of the projection of two members within 10". 
In any case, background contamination can only be esti- 
mated on a statistical basis and there is no single threshold 
that perfectly separate physical systems from projected 
ones. With this caveat in mind, we adopt the 2 a thresh- 
old in our analysis. 

In our survey, we therefore have identified 5 binary sys- 
tems with separation 0"8 < p < 10" out of 22 protostars 
in Taurus and 12 companions to 10 of our 41 sources in 
Ophiuchus. This corresponds to a companion star fraction 
of 23±9% in Taurus, 29±7% in Ophiuchus and 27±6% 
in the 1 10-1400 AU range for the combined sample. Our 
detection limit is not uniform across the separation range 
probed in this study and we may therefore be missing a 
few close, faint companions to our targets. We do not at- 
tempt to correct for this, as one would have to rely on 
assumptions regarding the flux ratio and separation dis- 
tributions of protostellar binary systems, which arc cur- 
rently unknown. Consequently, the number of companions 
we estimate might slightly underestimate the actual one 
over the separation range investigated here. In addition, 
we find one candidate companion in Taurus which is pos- 
sibly physically bound, though just at the 2(t level, to a 
YSO in the range 10-30". 

4. Discussion 

The survey we have conducted yields a statistically signif- 
icant estimate of the multiplicity frequency for a nearly 
complete population of embedded YSOs in Taurus and 
Ophiuchus. We discuss below the implications of this re- 
sult regarding the fragmentation of pre-stellar cores as 
they collapse ('5 14.1(1 . We then investigate whether prompt 
dynamical evolution of multiple systems occurs during the 
protostellar embedded stage, on a timescale of > 10^ yr, 
by comparing the statistics of multiplicity we derived for 
Class I and flat spectrum sources to that of even younger, 
more deeply embedded Class sources. We further com- 
pare the frequency of embedded multiple systems to that 
of Class II and III T Tauri stars in the same clouds, in 



order to investigate the dynamical evolution of multiple 
systems on a longer timescale of > 10^ yr fS I4.2() . Finally, 
we confront the predictions of various classes of models 
dealing with the formation and early evolution of multiple 
stellar systems to the statistics of multiplicity of embed- 
ded sources fS l4.3|l . 

4.1. Multiple protostellar systems and the 
fragmentation of collapsing cores 

The primary result of this study is that, over the separa- 
tion range 1 10-1400 AU, embedded YSOs present similar 
apparent multiplicity rates in the Taurus and Ophiuchus 
star-forming regions despite the known differences be- 
tween these two regions (pre-stellar core properties and 
stellar density). The size of prestellar cores being ^ 3 
times larger in Taurus than in Ophiuchus, the compari- 
son of the properties of binary systems spanning the same 
physical separation range may introduce a bias if core frag- 
mentation is a scale- free process (e.g., Sterzik, Durisen 
& Zinnecker 2003). We therefore also compare the 110- 
1400 AU multiplicity rate of Ophiuchus (29±7%) to that 
over the 330-4200 AU separation range in Taurus. The 
latter is 23±9% (5 companions out of 22 targets), sta- 
tistically indistinguishable from the Ophiuchus rate. This 
alternative comparison reinforces our conclusion that the 
frequency of companions among our two samples of em- 
bedded YSOs are identical in the 2 star-forming regions. 
This in turn suggests that, in spite of initial differences 
between Taurus and Ophiuchus pre-stellar cores, fragmen- 
tation appears to produce very similar properties for their 
multiple systems at the protostellar stage. 

Both samples of embedded YSOs we studied are nearly 
complete. Nevertheless, they exhibit a notable difference: 
roughly 25 % of the primaries in Ophiuchus have Lhoi < 
0.1 Lq, while only one such source is included in the 
Taurus sample. These low luminosity objects may either 
be extremely low mass objects (e.g. candidate embedded 
brown dwarfs) or embedded sources with reduced accre- 
tion rates. In any event, it might not be fair to compare 
the multiplicity properties of these two samples if they 
are made of intrinsically different populations (different 
mass function and/or evolutionary stages). The median 
luminosity of both samples is about 0.5 Lq. We compared 
the multiplicity rate of the faintest half of our samples to 
that of the brightest half. Overall, we find a marginally 
higher multiplicity rate among objects with Ltoi > 0.5 
(11 companions to 33 sources) than among less luminous 
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ones (6 out of 30) , with both molecular clouds suggesting 
the same trend although not at a statisticaly significant 
level due to the small sample sizes. Hence, irrespective 
of whether we include all embedded YSOs in our sample 
or not, the overall companion star fraction over the 110- 
1400 AU range is on order of 20-30 % in both star-forming 
regions. 

The observed proportion of companions to protostars 
we find is relatively high, especially if one considers that 
hardly more than an order of magnitude in separation 
is probed in this study. In the same separation range, 
the companion star fraction of field G dwarfs amount to 
only 14% (Duquennoy & Mayor 1991). This comparison 
is relevant if the distribution of separation for protostellar 
systems is as wide as that of mature systems. Kroupa & 
Burkert (2001) concluded from N-body simulations that 
stellar dynamical interactions in embedded clusters can- 
not significantly broaden the initial distribution of orbital 
periods, so that the wide range of separations observed 
for mature systems must result from the fragmentation 
process itself. At an earlier stage, accretion of material 
from the infalling envelope on a "seed" protobinary may 
also change the initial separation of the binary. Through 
SPH simulations. Bate (2000) has shown that this phe- 
nomenon is largely dependent on the (unknown) initial 
conditions. While we may expect accreting binary systems 
to get wider with time on average, we note that the excess 
of close companions in the separation range studied here 
is almost twice as large as the frequency of systems with 
larger separations for field G dwarfs. Hence, the excess 
companion star fraction we find over a restricted range of 
separation in embedded sources suggests that binary and 
multiple systems are indeed a very frequent outcome of 
the dynamical fragmentation of prestellar cores as they 
collapse. 

The similitudes of the multiple systems in the 2 molec- 
ular clouds extend beyond their frequency, and include as 
well their separation distribution and K band flux ratios, 
which both evenly sample the ranges probed in this sur- 
vey (Figure 13) • Protostellar multiple systems appear to 
span a larger range of flux ratios than the more evolved 
T Tauri (Class II) binaries, which are usually limited to 
AK < 3mag (e.g., Leinert et al. 1993). This may be be- 
cause the observed K band flux of an embedded YSO is 
mostly driven by its accretion rate and line of sight ex- 
tinction, while the central source's luminosity is usually a 
small fraction of the total system's luminosity (Bontemps 
et al. 2001). This is not the case for T Tauri stars, in which 
near infrared flux ratios more reliably trace mass ratios. 
Therefore, we cannot use the distribution of flux ratios 
of our sample of YSOs to infer the underlying mass ra- 
tio distribution for these objects, but can merely conclude 
that the relative accretion rates onto the components of 
the multiple protostellar systems appear similar in Taurus 
and in Ophiuchus, in spite of initially different prestellar 
core structures. 



4.2. The evolution of protostellar multiple systems 

The frequency and orbital properties of protostellar sys- 
tems is expected to evolve on various timescales due to a 
number of processes, such as dynamical decay of initially 
unstable configurations (~10'*~^yr, e.g., Delgado-Donate, 
Clarke & Bate 2003; Goodwin, Whitworth & Ward- 
Thompson 2004a) or disruptive encounters with other sys- 
tems if located in dense environments (~10^~^yr). In this 
context, a somewhat surprising result is the similar fre- 
quency and properties we find for protostellar systems in 
both Taurus and Ophiuchus. This suggests that their early 
dynamical evolution (<10^yr) has been remarkably simi- 
lar, in spite of both different initial pre-stellar core struc- 
tures and different environmental conditions, the stellar 
density in Ophiuchus being 2 to 3 orders of magnitude 
larger than in Taurus (e.g., Allen et al. 2002). 

Searching for possible evolutionary effects during the 
embedded phase, we first separate Class I from fiat spec- 
trum (FS) sources within our sample, as the latter are 
believed to be in an evolutionary state intermediate be- 
tween Class I and T Tauri objects, at least on a statistical 
basis. Combining the two star-forming regions in order 
to avoid small sample sizes, we find 9 companions among 
28 Class I sources and 8 companions among 35 fiat spec- 
trum primaries. This yields a companion star fraction of 
32±9% for Class I and of 23±7%, for FS sources over the 
110-1400AU separation range, a hardly significant evolu- 
tionary trend. 

Our source selection is, however, exclusively based on 
the near- to mid- infrared spectral index, which may not 
be the most appropriate tracer of a YSO's evolutionary 
status. For instance, T Tauri stars that are observed with 
an optically thick edge-on disk (such as IRAS 04158-1-2805 
in our sample) can have a spectral index typical of a Class I 
source. A more reliable indicator of the evolutionary sta- 
tus of an embedded YSO is the presence of an extended 
envelope, which characterizes a bona fide protostar. Such a 
protostellar envelope can be directly traced by its thermal 
millimeter continuum emission (e.g. Andre & Montmerle 
1994, Motte & Andre 2001) and often indirectly traced 
by an extended near-IR nebulosity (Park & Kenyon 2002 
and this study). The presence/absence of such an enve- 
lope is listed for all sources in our sample in Tables^and 
Twelve out of 22 YSOs in Taurus and 15 out of 41 
YSOs of Ophiuchus exhibit evidence for a protostellar en- 
velope. Combining the 2 regions, we find 12 companions 
to the 27 true protostars surveyed giving a companion 
star fraction of 44±10% in the 110-1400 AU range (see 
Table^J. In contrast, the more evolved sources in our sam- 
ple, which have already dissipated their envelope, have a 
much lower companion star fraction: 14±6 % (see Table0}. 
Hence, adopting the presence of an envelope as a qualita- 
tive diagnostic of youth, we find that young self-embedded 
sources tend to be more often in multiple systems than 
more evolved ones. Tfiis evolution of tiie multiplicity rate, 
seen liere at ttie 2.6a level, suggests a significant reduction 
of the number of wide multiple systems over a timescale of 
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Table 4. Observed companion star fraction over the separation range 110-1400 AU in Perseus, Taurus- Auriga and 
Ophiuchus star-forming regions (SFRs) for Class sources, Class I sources with (env-|-) and without (env-) an extended 
millimetric envelope, and for Class II-III T Tauri stars in the same clouds. For comparison, the multiplicity rate for 
G-type field dwarfs over the same separation range is 14% (Duquennoy & Mayor 1991). References: 1 - Looney et al. 
(2000); 2 - this study; 3 - Mathieu (1994); 4 - Barsony et al. (2003). 





SFR Nco„^p/Np„m CSF Ref. 

(%) 


CSFtot 
(%) 


CSFtot 
(%) 




Class 


Perseus 1 / H ot^+iQ i 
Ophiuchus 2 / 2 ^ 


25±13 


38±8 


Embedded sources 
with mm envelope 


Class I env-f 


Taurus 5/12 41±14 
Ophiuchus 7 / 15 47±13 


44±10 


Class I env- 


Taurus 0/10 
Ophiuchus 5/26 19±8 


14±6 


22±3 


Optical/IR sources 
without mm envelope 


Class II-III 


Taurus 24 / 100 24±4 3 
Ophiuchus 21 / 88 24±5 4 


24±3 



^ 10^ yr, which could conceivably result either from dis- 
ruptive encounters with other systems or from the internal 
decay of initially unstable systems. 

If this trend can be extrapolated back to even earlier 
stages of evolution, one may expect to find an even larger 
fraction of multiple systems among Class sources. The 
investigation of 14 YSOs driving large scale outflows led 
Reipurth (2000) to claim a multiplicity frequency in the 
range 79-86%, much higher than found here for Class I and 
flat spectrum sources. However, Reipurth's (2000) sample 
contains sources that span the full range from Class pro- 
tostars to T Tauri stars and some of the multiple systems 
have separations outside the range studied here. Moreover, 
if giant Herbig-Haro flows are actually triggered by the de- 
cay of unstable protostellar multiple systems, as Reipurth 
(2000) argues, this sample is obviously strongly biased to- 
wards the detection of such systems. Therefore, the de- 
rived multiplicity frequency of ~ 80% probably does not 
represent the true fraction of multiple systems in deeply 
embedded YSOs. A more meaningful comparison can be 
made with the results of the 2.7 mm continuum interfer- 
ometric survey of embedded YSOs performed by Looney, 
Mundy & Welch (2000) with a spatial resolution similar to 
ours. Their sample contains 12 "separate envelope" Class 
sources (see their Table 5, assuming that the members 
of the SVS13 group are all Class objects and classify- 
ing SSV13B as a separate envelope source), i.e., sources 
which would be considered as distinct primaries in our 
study. Among these, they find 3 companions in the separa- 
tion range 1 10-1400 AU corresponding to their definition 
of "common envelope systems" , to which their millimeter 
interferometric observations are mostly sensitive. In this 
range, the companion star fraction of their Class sam- 
ple is thus 3/12 = 25±13%, including 2 companions to 2 
sources in Ophiuchus. Such a fraction is consistent with 
the multiplicity rate we derive for Taurus Class I sources 
with extended envelopes (see Table 0J). Unfortunately, the 
poor statistics stemming from the small sizes of the sam- 
ples and the differing separation sensitivities of the mil- 
limeter interferometric survey and the present near-IR in- 
frared survey prevent us from drawing any strong conclu- 



sion about the evolution of the multiplicity rate between 
Class and Class I sources on a timescale <10^yr. 

Looking further downstream in the evolution of young 
multiple systems, we can also compare the multiplicity 
frequency we derived for embedded sources in Taurus 
and Ophiuchus to the fraction of multiple systems found 
among more evolved Class II and Class III T Tauri stars in 
the same regions. For the latter groups, we used the lists 
of multiple systems from Mathieu (1994) in Taurus and 
Barsony, Koresko & Matthews (2003) in Ophiuchus, which 
likely contain most O'.'8-IO" companions, as those can be 
found from direct infrared imaging. Beyond 1", those sur- 
veys are likely to be complete and they actually contain 
systems with flux ratios as large as 4 mag at least, thus 
similar to our detection limits. In the separation range 
considered here, there are 21 companions to 88 optically 
revealed T Tauri stars in Ophiuchus and 24 companions 
to 100 primaries in Taurus (see Figure^. The two pop- 
ulations have a companion star fraction of 24±3%, i.e., a 
factor of about 2 lower than that of Class I sources with 
envelopes (44±10%). The much lower multiplicity frac- 
tion of Class II-III sources in the 2 clouds compared to 
that of the youngest Class I sources surrounded by mil- 
limeter envelopes (a 1.9(t difference) seems to confirm the 
similar trend reported above between young and evolved 
embedded Class I sources. 

Table^lsummarizes the statistics of multiplicity among 
embedded sources and young PMS stars in Taurus and 
Ophiuchus. The most significant result, yet still rely- 
ing on relatively small samples, is provided by the com- 
parison between deeply embedded Class and Class I 
sources surrounded by extended millimeter envelopes on 
the one hand, and more evolved, envelope-less Class I-IT 
III sources, on the other. The companion star fraction is 
seen to decrease from 38±8 % in the youngest embedded 
sources (age< 10^ yr) down to 22±3 % in pre-main se- 
quence objects (age ^10^ yr). Overall, we thus find ev- 
idence for an evolution of the multiplicity fraction dur- 
ing the embedded phase, with a decrease of multiple sys- 
tems over a timescale of 10^ yr at most, significant at the 
2(7 level. The comparison between Class I sources with 
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and without extended envelope suggests that the actual 
timescale is on the order of ~ 10^ yr. 

We caution, however, that this result is limited by rel- 
atively small number statistics and strongly depends on 
the reliability of the evolutionary classification adopted 
for the YSOs in our sample. While the mass and/or opti- 
cal depth of the envelope is generally recognized as being 
a robust tracer of the circumstellar evolution of embedded 
YSOs (e.g. Myers et al. 1998; Andre, Ward-Thompson, & 
Barsony 2000), it is more difficult to obtain reliable es- 
timates of envelope properties in cluster-forming regions 
such as the Ophiuchus main cloud than in regions of 'iso- 
lated' star formation such as Taurus. In Taurus, the rel- 
atively large spatial extent of prestellar cores and pro- 
tostellar envelopes (~ 18 000 AU) makes it easy to dis- 
tinguish between bona-fide protostars with envelopes and 
more evolved Class I sources (Motte & Andre 2001). In 
Ophiuchus, however, the smaller prestellar fragmentation 
lengthscale leads to more compact protostellar envelopes 
which are more difhcult to spatially resolve with single- 
dish millimeter telescopes. The current classification of 
Ophiuchus Class I sources in objects with or without ex- 
tended envelopes (see Table|3 Andre & Montmerle 1994; 
Motte et al. 1998) is thus more uncertain than in Taurus, 
pending the results of interferometric measurements in 
progress. If we consider the Taurus cloud alone, the de- 
crease in the multiplicity fraction from self-embedded 
protostars to envelope-less objects becomes apparently 
stronger, although the sample size is then only barely sta- 
tistically significant (see Table 0)). In this framework, our 
results suggest a steep decrease of the fraction of multi- 
ple systems occurring during the protostellar embedded 
phase as the envelope is being dissipated on a timescale 
of-lO^yr. 

4.3. Confrontation with model predictions 

The above results indicate that the fraction of wide (110- 
1400 AU) protostellar multiple systems is initially high 
and seems to decrease on a timescale of ~ 10^ yr during the 
Class I phase in the Taurus and Ophiuchus clouds. How do 
these results compare with models of the formation and 
evolution of multiple systems ? 

Assuming that the distribution of orbital periods of 
embedded multiple systems is at least roughly as wide 
as that of T Tauri and field binaries (Kroupa & Burkert 
2001), the multiplicity fraction we derive over a restricted 
range of separation suggests that, on average, each embed- 
ded YSO has at least one companion. This result is consis- 
tent with a class of models which rely on the assumption 
that all stars are initially in binary systems and investigate 
their subsequent dynamical evolution (e.g., Larson 1972; 
Kroupa 1995b). Such models predict that the initially high 
fraction of multiple protostellar systems significantly de- 
creases on a timescale of 1 Myr or less in dense protostellar 
clusters, such as ONC, due to frequent disruptive encoun- 
ters between systems (Kroupa, Aarseth & Hurley 2001), 



whereas the dynamical evolution is much milder in loose 
star-forming regions, such as Taurus, where gravitational 
encounters are rarer. 

The high fraction of embedded multiple systems we 
report here for Taurus is consistent with the predictions 
of these models for low density regions. The small amount 
of dynamical evolution expected to occur in Taurus will 
take place very early on (< 10^ yr), when the initial proto- 
stellar aggregates still have a density of order of 10^ stars 
pc~^ (Kroupa & Bouvier 2003). Disruptive encounters be- 
tween protostellar systems at this stage may result in the 
decreasing multiplicity fraction we observe in Taurus dur- 
ing the Class I phase on a timescale of ~ 10^ yr (Table 0)). 
Past this early stage, the stellar aggregates quickly ex- 
pand on a timescale of ~ 10^ yr to reach the currently 
observed densities of order of 10-100 stars pc~^ and the 
Taurus population is not expected to dynamically evolve 
any more. Thus, the higher multiplicity rate of both em- 
bedded YSOs and T Tauri stars in Taurus with respect 
to field stars suggests that loose star-forming regions are 
not the dominant mode of star formation, i.e., that only 
a small proportion of all field stars form in such areas. 

On the other hand, Ophiuchus is a young embed- 
ded cluster similar to the "dominant-mode cluster" of 
Kroupa's (1995a) models, from which most field stars are 
expected to arise (Clarke, Bonnell & Hillenbrand 2000; 
Carpenter 2000; Adams & Myers 2001). In Ophiuchus, 
as for Taurus, the initial decrease of multiplicity fraction 
we observe during the Class I embedded phase may be 
accounted for by disruptive encounters occurring in high 
density protostellar subclusters on a timescale of ^lO^yr 
(Bonnell, Bate & Vine 2003). Unlike Taurus, however, the 
remaining excess of binaries observed among Class II-III 
T Tauri stars in Ophiuchus at an age of ~10^yr has still 
to be reduced by a factor of 2 to be consistent with the 
frequency of field binaries. The current average density of 
about 10^ stars pc""^ in the Ophiuchus cloud translates 
into a crossing time of order of 1-2 Myr (e.g., Kroupa 
1995a), which may indicate that the T Tauri population 
is not fully dynamically evolved yet in this cloud. 

Together with observations of embedded multiple 
YSOs in an even denser cluster such as Orion, N-body 
simulations of the evolution of an Ophiuchus-likc cluster 
would provide a crucial test on the assumption that core 
fragmentation produces a universal set of initial proper- 
ties for protostellar multiple systems which subsequently 
evolve depending only upon environmental conditions. We 
note however that the similar frequency, properties and 
evolution we derive for embedded multiple systems in both 
Taurus and Ophiuchus seems to argue against large-scale 
environmental conditions playing a crucial role in shaping 
the properties of protostellar systems. Instead, the similar 
evolution of protostellar multiple systems in star-forming 
regions whose average properties (e.g., mean stellar den- 
sity) are so different may point to local conditions pre- 
vailing in e.g. small-scale embedded subclusters as being 
determinant (Bonnell et al. 2003; Sterzik et al. 2003). 
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The observed high frequency of multiple systems 
among embedded YSOs indicates that most prestellar 
cores undergo fragmentation into at least two objects dur- 
ing their collapse. It is often assumed that prestellar cores 
actually experience multiple fragmentation as they col- 
lapse thus leading to the formation of several protostellar 
seeds, i.e., small-N protostellar aggregates with TV = 3-10, 
typically, on a scale of a few 100 AU (e.g. Burkert, Bate 
& Bodenheimcr 1997; Stcrzik & Durisen 1998; Dclgado- 
Donate et al. 2003; Goodwin et al. 2004a). Global simu- 
lations of cloud collapse support this scenario and predict 
the formation of even higher density subclustcrs consisting 
of 10-50 protostellar seeds (Klessen, Burkert & Bate 1998; 
Bate, Bonnell & Bromm 2003; Bonnell et al. 2003). Small- 
N protostellar aggregates as well as high density subclus- 
tcrs undergo rapid dynamical evolution on a timescale of 
lO^-s yr (Sterzik & Durisen 1998; Delgado-Donate et al. 
2003). The dynamical decay of initially unstable groups is 
thus expected to eject most of the lowest mass protostellar 
seeds, with a typical ejection velocity of < 1 km/s, while 
those remaining in the cloud core form stable binary and 
higher order multiple systems. 

At an age of about 10^ yr, we find no evidence that em- 
bedded YSOs in Taurus and Ophiuchus belong to small- 
N clusters or high density subclusters on a scale of a few 
100 AU. No triple embedded YSO is found in Taurus, and 
only 2 in Ophiuchus, compared to a total of 13 binary 
YSOs detected over the separation range 1 10-1400 AU. 
It may be argued that dynamical evolution of small-N 
clusters has already taken place by that time, ejecting 
most of the initial protostellar seeds at larger distances 
(^1000 AU) and leaving only the resulting stable sys- 
tems. However, Allen et al. (2002) investigated the spatial 
distribution of prestellar cores and YSOs in Ophiuchus 
and concluded that the YSOs had neither concentrated 
nor dispersed significantly since their formation in these 
cloud cores. Numerical simulations of the dynamical decay 
of small-N clusters including accretion onto the protostel- 
lar seeds and circumstcUar disks (e.g. Delgado-Donate et 
al. 2003, Goodwin et al. 2004a) predict the formation of an 
equal number of tight binaries and higher order (triple and 
quadruple) systems. This prediction also seems to contrast 
with our results, although it may simply indicate that a 
fraction of our sample primaries actually arc tight bina- 
ries with a separation less than 100 AU and have therefore 
gone unresolved in this study. Pending a higher resolu- 
tion study of these embedded sources (Duchene et al., in 
prep.), a comparison of the hierarchy of multiplicity as 
derived from observations to model predictions remains 
somewhat inconclusive. 

For the time being, on the sole basis of the high fre- 
quency of wide binaries among embedded Class I and flat 
spectrum sources, we cannot firmly distinguish between 
the formation of unstable small-N aggregates and a core 
fragmentation process that only creates 2 or 3 seeds, 
i.e. stable multiple systems without further internal dy- 
namical evolution. In any event, the similar multiplicity 
rate we derive for embedded YSOs in both Taurus and 



Ophiuchus in spite of initially different prestellar core sizes 
and widely different average stellar density between the 
2 star- forming regions indicates that, regardless of the 
formation channel, the outcome is fairly insensitive to 
both initial and environmental conditions in this range. 
Whether this is accounted for by current models remains 
to be seen (e.g., Delgado-Donate, Clarke & Bate 2004; 
Goodwin, Whitworth, & Ward-Thompson 2004b; Sterzik 
et al. 2003). 

5. Conclusions 

We have conducted a survey of visual companions to 63 
embedded YSOs in the Taurus-Auriga and Ophiuchus 
star-forming regions. The main results are as follows: 

— The companion star fraction of Class I and flat spec- 
trum embedded sources over a separation range 110- 
1400 AU is 23±9% in Taurus, 29±7% in Ophiuchus 
and 27±6 % for the combined sample. This fraction is 
twice as large as the freqiiency of companions to G 
dwarfs in the solar neighborhood, in the same separa- 
tion range. Assuming that the distribution of orbital 
periods of protostellar multiple systems is about as 
wide as that of mature systems, this high multiplicity 
fraction suggests that each embedded YSO has, on av- 
erage, at least one companion. The fragmentation of 
prestellar cores as they collapse thus appears to be a 
frequent occurrence indeed. 

— In spite of the known differences in the properties 
of prestellar cores in Taurus and Ophiuchus and the 
widely different stellar density in the 2 star-forming 
regions, we find identical properties for their embed- 
ded multiple systems (frequency, separations, luminos- 
ity ratios). This suggests that the outcome of the frag- 
mentation of prestellar cores is fairly insensitive to 
both initial and environmental conditions, at least for 
the range of physical conditions encountered in these 
2 clouds. 

— Combining the Taurus and Ophiuchus clouds, the frac- 
tion of multiple systems appears to decrease on a 
timescale of order of ^ 10^ yr during the embedded 
phase, with the youngest embedded sources. Class 
and Class I YSOs still surromidcd by extended milli- 
metric envelopes, exhibiting nearly twice as many miil- 
tiple systems as more evolved Class I YSOs devoid of 
mm envelopes and Class II-III T Tauri stars. This sug- 
gests that rapid disruption of wide protostellar multi- 
ple systems occurs during the embedded phase, result- 
ing either from the internal decay of initially unstable 
systems or from disruptive encounters between stable 
systems in high density protostellar aggregates. 

— Over a scale of a few 100 AU, embedded YSOs at an 
age ~ 10^ yr in Taurus and Ophiuchus clouds are found 
to be in relative isolation and not assembled in subclus- 
ters. This seems to contrast with the assumption that 
protostellar collapse leads to small-N clusters (N=3- 
10) as well as with global simulations of cloud collapse 
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which predict the formation of even higher density 
subclusters on these scales. Evidence for such cluster- 
ing properties may have to be searched for on cither 
smaller scales (tight binaries) or larger ones (ejected 
protostellar seeds) if dynamical evolution has already 
occurred. 

Taurus and Ophiuchus clouds probe different modes 
of star formation, the former being a low density associ- 
ation harboring only low mass objects, while the latter is 
a dense cluster including at least one massive B star as 
well as a large population of low mass T Tauri stars. The 
results reported here unexpectedly suggest that the prop- 
erties and early evolution of embedded multiple systems 
in these two clouds are quite similar and therefore do not 
seem to depend much on initial and/or global environ- 
mental conditions. This may indicate that the properties 
of multiple systems are ultimately driven by dynamical 
processes acting on small scales and are thus shaped by 
local (rather than global) conditions which might not be 
very different between clusters and associations during the 
early stages of protostellar evolution. 
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Appendix A: TTauri stars in Ophiuchus 

In the course of our survey, we observed many objects 
that are members of the Ophiuchus star-forming regions 
but are in the more evolved classical and weak-line T Tauri 
(Class II and III respectively) stages. For completeness, we 
list here all the objects we observed (Table IA.1|) as well 
as the candidate companions detected within 10" of them 
(Table rO|l . Among the companions listed in Table FOl 
only 3 were previously known, namely the ISO 60-61 pair 
(listed in Bontemps et al. 2001) and the companions to 
ISO 95 and ISO 204, listed in Haisch et al. (2002). 

It can be noted that we detected a much larger fraction 
of candidate companions among this sample of T Tauri 
stars than among the protostars that form the primary 
sample of this study (51±6% vs. 30±6%). This is most 
likely a consequence of the fact that the protostars are lo- 
cated deeper in the molecular cloud and therefore they are 
surrounded by comparatively darker areas in our images. 
On the other hand, most T Tauri stars are located in ar- 
eas that probably contain many background stars that are 
not completely extincted, which leads to a larger number 
of projected visual pairs. This is confirmed by the gener- 
ally large flux ratios between the T Tauri stars and their 
apparent companions: almost 75% of all pairs listed in 
Table rO have AK > 5 mag while only 3 such pairs are 
found in Table |31 



16 



Duchene et al.: Multiple protostellar systems 



Table A.l. List of all T Tauri stars (Class II and III objects) in Ophiuchus that were detected in the vicinity of our 
target protostars. Objects whose name starts with "B" are from Barsony et al. (1997). The classification is determined 
from the infrared spectral indices following Bontemps et al. (2001). 
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Table A. 2. List of all companions to T Tauri stars in Ophiuchus that were detected during the course of our survey. 
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Notes - a: The companion to this source is ISO 61 (=GY103), a Class III object, 
b: The primary is saturated at K band; no observation is available in H band, 
c: The primary is slightly saturated. 
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